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Abstract. We review a recent theoretical development based on non-equilibrium Green's 
function formalism to study heat transport in nanomechanical devices modeled by phononic 
systems of coupled quantum oscillators driven by ac forces and connected to phononic reservoirs. 
We present the relevant equations to calculate the heat currents flowing along different regions 
of the setup, as well as the power developed by the time-dependent forces. We also present 
different strategies to evaluate the Green's functions exactly or approximately within the weak 
driving regime. We finally discuss the different mechanisms in which the ac driving forces deliver 
the energy. We show that, besides generating heat, the forces may operate exchanging energy 
as a quantum engine. 



1. Introduction 

The Green's function formalism provided the theoretical framework for many studies of 
electron transport in mesoscopic devices. Since the pioneer work by Caroli and coworkers, 
PQ the celebrated generalized Landauer-Biittiker formula by Meir and Wingreen [2] and the 
generalization to setups driven by the time-dependent voltages and fields of Refs. [21 H], several 
interesting applications can be found in the literature. Details of these historical developments 
and some important aspects of the subsequent progress have been reviewed in papers of this 
conference series [5j [6] . 

In the last years, transport experiments carried out in the presence of ac driving fields and 
the study of the related pumping phenomena [Zl El El QUI [HI EE2 Q3] motivated the development 
of strategies to formulate and efficiently solve the corresponding Dyson's equations in theoretical 
models for these setups. A rigorous formulation was presented in Refs. [13] for the problem of 
a mesoscopic ring threaded by a time-dependent magnetic flux. Such treatment was afterwards 
adapted to study particle transport in models of quantum pumps consisting in quantum dots 
where local ac voltages are applied, with [J5] and without many-body interactions [IB]. In the 
latter case, an explicit relation between the Green's function and the Floquet scattering matrix 
was presented in Ref. [17] . 

Microscopic models and the framework provided by non-equilibrium Green's functions also 
offer the adequate basis for the study of mechanisms for heat generation, energy transport 
and quantum refrigeration by electrons in mesoscopic structures under ac driving, [18] as well 
as the possibility of characterizing these processes by defining effective temperatures [19j . In 
the context of phononic systems, J-S Wang and coworkers proposed recently in Ref. [20] a 



treatment analogous to that of Ref. pQ to describe thermal transport induced by a temperature 
gradient under stationary conditions. The present contribution is devoted to review a recent 
related development, which focuses on the study of heat transport by phonons under ac driving 
|21| . One of the motivations is to investigate cooling mechanisms in these systems. This is 
because refrigeration by means of manipulating phonons, rather than electrons, has the appealing 
feature of being suitable for insulating as well as metallic systems. On the other hand, cooling 
nanomechanical systems is a very active field of research [22J 123 (Ml E51 [23 E3 [25] and it is 
interesting to explore alternative mechanisms. 

In a recent work, [25] it was shown that when a barrier (or perturbation) moves along a cavity 
connecting two gases of acoustic phonons, a sizable amount of heat can be transferred from the 
colder to the hotter gas. When the process is repeated as a cycle, the result is refrigeration until 
the colder reservoir reaches a minimum temperature determined by geometrical parameters of 
the device and the velocity at which the barrier travels through the cavity. That mechanism 
was afterward analyzed in the framework of a microscopic model which was solved by means 
of a nonequilibrium Green's functions treatment analogous to the one proposed to study heat 
transport in ac driven electron systems. In the subsequent sections we briefly review that 
theoretical treatment. We also present some approximations to study the problem in the limit 
of weak driving. Finally, we show that, besides refrigeration, another interesting mechanism 
suggested to take place in quantum pumps [18] and quantum capacitors [30], can be also realized 
in phononic systems. Namely, the exchange of energy between the driving fields. 

2. Background 

2.1. Microscopic Model 

The standard microscopic model for acoustic phonons consists of masses coupled by springs. 
The one-dimensional (Id) version of this model became popular in the context of heat transport 
since the pioneer work by Fermi, Pasta, and Ulam [31J and several subsequent studies devoted to 
understand the violation of the Fourier law in these systems [20J [321 [331 EH ESI EU] . Fourier law 
implies a linear drop of the local temperature along a system connected between reservoirs at 
different temperatures. It is the counterpart of Ohm law in the context of heat transport, and it 
is nowadays rather natural to accept that it should not be expected to hold through a finite chain 
of harmonic oscillators, along which heat propagates coherently. In fact, a similar effect is known 
to take place in the behavior of the local voltage along a mesoscopic structure placed between two 
reservoirs at different chemical potentials |37 [ 1381 [39]- I n these systems, the voltage drop takes 
place at the contacts and the concept of "contact resistance" was introduced to characterize this 
feature, while it vanishes along the structure through which electrons propagate coherently. It is 
remarkable that the studies on the origin of the lack of the Fourier law in phononic systems and 
the introduction of the concept of the contact resistance in electron systems followed parallel 
routes and the connection between these two phenomena has not been mentioned even in review 
articles [351 136] , More recently, heat pumping has been analyzed through a two-level system 
driven by an ac field coupled to phononic baths at different temperatures, [H] between semi- 
infinite harmonic chains with a time-dependent coupling, [42] and in anharmonic molecular 
junctions in contact to phononic baths |43| . 

We consider the one-dimensional (Id) system sketched in Fig. [H with a central finite chain 
of N molecules with identical masses m c coupled by springs of constant k^. This system is 
connected, at certain sites l a , to semi-infinite chains of masses m a coupled by springs of constant 
k a . These chains play the role of reservoirs, which we assume being kept at temperatures T a . 
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Figure 1. Sketch of the setup containing a central chain of acoustic phonons represented by a 
Id lattice of atoms with equal masses connected by springs. This system is in contact to several 
reservoirs at different temperatures and affected by driving forces acting locally at the different 
positions of the central chain. 

The first term H a represents the reservoirs; the second H c (t) describes the central chain. We 
assume that external time-dependent perturbations are applied at this system. The Hamiltonian 
reads 

N „2 AT-1 , N ,//.s 

H c (t) = Ef^ + E - * c ,m) 2 + E -T <u (2) 

1=1 Zm c 1=1 £ 1=1 £ 

where the last term of this expression describes the time dependent perturbation, which is 
represented by time-dependent elastic forces at different positions of the central chain. It is 
useful to define two characteristic frequencies in the problem: the frequency oj c = y/ko/m c , and 
the driving frequency f2o-(We shall work on units such that h = ks = 1). The contact between 
the central chain and the reservoirs is described by the Hamiltonian 

L ^ 

#cont = E ~§-( X <*,l ~ X c,lS- ( 3 ) 
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Notice that the site j = 1 of each reservoir couples to the site l a of the central region. The 



reservoir Hamiltonians are given by 



j=l L 



Pa,j ka , ^ .2 
2m a + 2 {Xa ' j Xa ' j+1) 



(4) 



In the limit of semi-infinite chains, N a — > oo, it is convenient to express the degrees of freedom 
of the reservoirs in terms of normal modes. For open boundary conditions, this corresponds to 
performing the following transformation, 
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The corresponding Hamiltonians transform into 



n = 0,...,N a . 
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2.2. Energy balance 

Following a procedure similar to that used in Ref . |18] , we consider the evolution in time of the 
energy density stored at a given elementary volume and write the corresponding conservation 
equation. Since we are treating a Id lattice, and our Hamiltonians involve only nearest-neighbors 
terms, the minimum volume that we consider is one that encloses nearest-neighbors sites, like 
the left box of Fig. If we denote by dE^i + i/dt the total energy flow exiting and entering a 
volume of the central chain without connection to reservoirs, the equation for the conservation 
of the energy for reads: 

= j? +1 (t) - J?® + m) + p l+1 (t), (9) 

where J®(t) denotes the incoming energy current from site / — 1 towards the site / and J^ 1 (t) 
the outgoing flow from site I + 1 to / + 2. We denote with a positive sign the flows pointing from 
left to right. These quantities are evaluated from 

dEii +x d[Hi + Hi + i}(t)) 

dt = -i{[H, H h i + H l+ljl+1 + H hl+1 \) + — , (10) 

where H\ ii are the matrix elements of H c (t) on the sites The explicit calculation cast 

J? = ^c,lX c ,l-l), (11) 
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for the current and 



PM = ^(4,i), (12) 



for the power developed by the force applied at the site I. 




Figure 2. Sketch of the setup indicating the construction employed to derive the expressions 
for the heat currents as well as the poser developed by the external forces. 



2.2.1. Heat flow into reservoirs We now focus on the right box of the sketch of Fig. [2j The 
variation in time of the energy stored in the reservoir a is 

^> = -i([H,H a ])=J$(t), (13) 

where we assume that no driving force is acting at the contact site of the reservoir. The current 
Ja(t) flows along the contact between the site l a of the central chain and the reservoir a and it 
is assumed to be positive when it enters the reservoir. Since we are dealing with phonons, this 
energy flow is a pure heat flow and its explicit expression reads 



■#(*) =E^Jjr^J S ^n)(^,n P c,l a )- (14) 



n 



Conservation of energy implies that the total average of the power invested by the external fields 
is dissipated into the reservoirs at a rate 

L N 

E 7 ? = £^ ( 15 ) 

a=l 1=1 

where 

P t = - f dtPi(t), (16) 

T JO 



with r = 27r/f2oj is the power developed by the external fields at the site I averaged over one 
period, while 

jj = l [ T dtjg(t), (17) 

T JO 

is the dc component of the heat current at the reservoir a. Notice that a net amount of work 
has to be invested in order to pump heat from one reservoir to the other. In the next sections we 
show that carrying out an analysis similar to that developed in Ref. [18] leads to the conclusion 
that the rate at which the total work done by the external fields is dissipated as heat flowing 
into the reservoirs is proportional to $7q. 



3. Non-equilibrium Green's function approach 

3.1. Relevant Green's functions 

The theory of non-equilibrium Green's functions, which was developed independently by L. P. 
Kadanoff and G. Baym |44| . Schwinger |45j and Keldysh [26] has been employed many times 
in the context of electron transport, as mentioned in the introduction. Details can be found in 
that literature. In the present case, the general procedure is completely analogous. The relevant 
Green's functions are the bigger and lesser functions 



£>>(1,1') 
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where X(l) = y/fh~ixi(ti) is a phononic operator and 1, 1' is a schematic notation that labels 
the spacial and time coordinates. The corresponding retarded and advanced Green's functions 
read 



D R {l,l') 
D A (1A') 



8(t 1 -t 1 ,)[D>(l,l')-D<(l,l')} 

e(t 1 ,-t 1 )[D<(i,i')-D>(i,i% 



where the four Green functions satisfy the following relations D A (1,1') 
£><(!, l') = -[£>>(!', 1)]* and D R (l, 1') - D A (1, 1') = D>(1, 1') - D<(1, 1'). 
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3.2. Dyson's equation 

We now turn to derive the Dyson's equations to evaluate the Green's functions in our problem. 
The operators x and p satisfy canonical commutation relations: 

[x u x v ] = 0, [xi,pi>] = iSijr. (21) 

For coordinates 1,1' within the central chain, the lesser and greater Green's functions read 

D<,(t,t') = im c {x c j,(t')x c> t(t),) (22) 

-^'(M') = im c (x c> i(t)x c ^(t r )}, (23) 

while the retarded Green's function is 

£>5'(M') = -ie{t-t')m c ([x c ,i{t),x c , v {t')]). (24) 

Now we proceed to the calculation of the equation of motion of the retarded Green's function. 
To this end we must take the second derivative respect to t' , and use the Eherenfest's theorem. 
The result is 
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(25) 
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where l a refers to the site of the central chain connected to reservoir a and a\ the first site of 
the same reservoir. In turn, the integral representation for the Dyson's equation corresponding 
to the function D R reads 
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with 77 = + stands for the equilibrium Green's function within the uncoupled reservoir a. If we 
replace the Eq . (I27P into Eq. ([25j) . we can define the self energy as 
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- t') + D R (t, t')M{t') - I dtib R (t, h)t R (h, t') = s(t - 1') 

with M(t) = M(°) + MW(i), where 



(a more explicit expression is presented in Appendix A), and Eq. (|25p can be written in the 
simpler matrix form 

(30) 
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and Eg,(t,f) = E«=l V^aCM')- 

On the another hand, the Dyson's equation for the lesser function is obtained by recourse to 
Langreth's rules 



£>M'(M') = Y,J dt iJ d ^ D L^ t i)^(ti-t 2 )D^ ll (t 2 ,t'), 
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and 
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where Eq. (|33p corresponds to coordinates along the central chain while Eq. (|34p corresponds to 
one coordinate along the chain and the other on the reservoir. The advanced Green's function 
of the uncoupled reservoir is d^ a (ti,t') = [d R a (t',ti)]*, while the corresponding lesser function 



(34) 



reads 



(35) 



with n a (w) = \ j{& 1 l T < x — 1) being the Bose-Einstein distribution function, which depends on the 

] the 

(36) 



temperature T a of the reservoir a. In terms of the spectral function defined in Appendix A the 
Fourier transform of the lesser Green's function reads 



S<(w) = m tt (w)r n (u). 



3.3. Solution of Dyson's equation 

To obtain the retarded function of the central region, D R , one has to solve the differential 
equations (|30p. For the case of terms that depend harmonically on time, it is convenient to 
follow the strategy proposed in Refs. |14[ [T6] for electron systems and recently adapted to 
phonons in Ref. [21 J. We summarize the main steps in the present section. It is convenient 
to perform the Fourier transform with respect to the "delayed time" in the retarded Green's 
function 

D R (t,u)= [ + °° dt'e^+^-^M^t'), (37) 



and to represent the time-dependent perturbation of the Hamiltonian in terms of its Fourier 
expansion 

K 

M«(t) = M^e - **" *. (38) 

k=-K,k^0 

Substituting in the Dyson's equation (f3"U|) results 

+oo 

D R (t,u) = D ( -°\u)+ e-^^it^ + k^M^D^ico), (39) 

kj^O,— oo 

where 

L> (0) (w) = [io 2 I - M (0) - E*(w)l ^ , (40) 

corresponds to the stationary component of the retarded Green's function of the central chain 
connected to the reservoirs but without the effect of the time-dependent perturbations. Notice 
that because of the periodic structure of the time dependent part of the Hamiltonian, the 
retarded Green's function is also periodic in time with period r. For this reason it is useful and 
natural to represent this function in terms of a Fourier series as follows 

+oo 

D R (t,u)= e- iknot V(k,u) (41) 

k=— oo 

where the functions T>(k, uj) are also known as Floquet components. The exact solution for the 
set of coupled linear equations Eq. ([39j) can be obtained numerically by following the procedure 
introduced in Ref. [47] for fermionic systems driven by ac potentials. Alternatively, it is 



possible to solve these equations approximately in the limits of weak amplitudes of the 

time-dependent perturbations and in the limit of low frequencies Qq. We briefly present these 
approximate treatments in the next subsections. 



3.4- Perturbative solution of Dyson's equation 

In most cases, the set of linear equations of Eq. (|39p must be solved numerically. For this reason, 
it is convenient to carry out a systematic expansion in powers of M^ l \ to obtain analytical 
expressions. When these amplitudes are small (compared with the energies of the Hamiltonian 
independent of time), a perturbative solution in these parameters is a good approximation. 
Assuming that the external perturbations M^\t) contain K harmonics, as expressed in Eq. 
(|38p . the Green's function evaluated up to the 1st order in the amplitudes Mu^ is 



K 

D R (t,cu)~ e- ikQot V(k,oj), 



(42) 



with 

P(0,oj) = P(w), (43) 
V(±k,u) = D ( - 0) {oj±kn )M ( ± ) k DW{uj), k = l,...,K. (44) 

3.5. Low driving frequency solution of Dyson's equation 

Let us now consider the Dyson equation in Eq. (|39p in the limit of low driving frequency CIq. 

1st order A solution exact up to 0(Qo) can be obtained by expanding Eq. (f3"9"j) as follows: 

D(t,u) ~ D^(u) +D(t,u)M ( - 1 \t)D ( -°\uj) + 

(45) 

We define the frozen Green function 

D f (t,u)= p^-Pd" 1 , (46) 
in terms of which the exact solution of the Dyson equation at (D(CIq) reads 

D^(t,uj) = D f (t,u) + id w D f (t,u) dM ^® £><®(u). (47) 

2nd order To obtain the solution exact up to 0(Qq), we consider the following expansion of 
equation Eq. (|39|): 

D(i,cj) ~ b^\u) + b{t,u)M^{t)D^{u) 

at 

4^( ( ,.)«V>M. (48, 

The solution exact up to 0(fij§) is 

£( 2 )(t, W ) = D f (u) + id w D( 1 \t,u) dM ™® D f (u) 

~\dlDf{tM^^H-)- (49) 

One can obtain this Green function numerically by first discretizing the time in the interval 
< t < t, then solving Eq. (|49p for each one time, and finally evaluating the Fourier transform 
to obtain the Floquet representation with which the dc component of the heat current can be 
calculated from Eq. (|55p. 

4. dc heat current 

In this section, we present the steps to calculate the net heat current flowing into a given 
reservoir. We start by writing the time-dependent heat current given in Eq. (|14p as follows 



Jg(t) = £^viv r TT 8in( *" ) t5?^ <x °'" (i ' )x, " (<)) (50) 



We now substitute Eq. flMJ) and the representation of Eq.(J3T]) to write the dc component of the 
above current as 



L +00 

E E ^ + kn )\v lath (k,u)\ 2 i:<(u)d^(u J + kn )} (52) 
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where k n ^ a = ( m ^J ]^~+T sin2 (<7n)- Finally we obtain 
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(53) 

An alternative representation to this current can be obtained after recasting the first term of 
the above equation by recourse to the identity 

V hv {k,uj)-V* Vjl {-k,uj + kn Q ) = -iY^Y. V hi a {k + k\u-k'Sl Q ) 

k' a 

xr a (cj- k'noMjtk',^- k'n ). (54) 



This leads to the following equation for the dc heat current, 
— n L + °° r+°° duo 

Ja=Y, E J_ ^ + kn )[n p (Lj)- n a (cj + kn o )]T a (oj + kn o )T (oj)\V lajl(! (k,oj)\ 2 . (55) 
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This equation constitutes a generalization to the case of ac driving of a Landauer-Biittiker 
formula for the heat transport in phononic systems (see Ref. |20j) 



/ + OO 
7-w^T^(w)[^(w) -n a (u)], (56) 
-°° Z7r /3=1 

being 7^/3 (w) = \D(^ t (oj)\ 2 T a (uj)Tj3(oj), the transmission function between the reservoirs a 
and j3. Equation (|55|) indicates that a net heat current may exist even in the absence of a 
temperature difference between the reservoirs. Such a flow is a consequence of the pumping 
with the ac forces and contains, in general, an incoming component which accounts for the work 
done by the ac fields and which is dissipated into the reservoirs. As shown in Ref. [21] under 
certain conditions, a net current between the reservoirs can be established. This current can go 
against a temperature gradient, therefore allowing for cooling. 



5. Mean power developed by the external forces 

We now turn to analyze another interesting aspect of the heat transport, which are the 
mechanisms used by the driving forces to transfer a net amount of energy to the phononic 
system. For sake of simplicity we restrict our analysis to the case of driving forces that oscillate 
with a single harmonic component and a phase-lag, of the form 



k[{t) = K cos(n t + 5i). 



(57) 



The mean power developed by the force acting on the site I corresponds to calculating the 
dc component of the time-dependent power of Eq. (|12p . In terms of Green's functions can be 
written as follows 

-i r*a*K«)n<. 



Pl = / dt ^l D <(t,t). (58) 

rm c Jo at 

According to the Dyson's equation (|33p and in terms of the Floquet-Fourier representation of 
the Green's function ()4ip the above equation reads 



- = _Qo ^ f^^M^h'V^k^MVl^k + k'^), (59) 



a k,k'' 

For the case of ac forces of the form (|5Tj) this equation results 



p i = ^rEE / "^MM^M^Ii-i,^ (60) 



k a=\ ' 



-V lu {k^)T a {oj)Vl loi {k + l^)e-^}. (61) 

A simple analytical result can be derived for weak driving, in which case we can use the 
perturbation treatment of Sec. 13,41 This leads to a result for the mean power which is exact up 
to 0{Kq) and reads 

_ 9nK 2 N r +0 ° 

P l ~ TT^ — E / du}lm i l COS ( 5 l ~ 5 i) ~ i Sin (^ " 5 3j\ ( 62 ) 

%j (u;)D^(u - 9 ) + 7 fc(a,)Z?£>(o; + 9 )} } , (63) 



with 



7!» = E n a ( W )L»£*(a;)r a (a;)^; ) a (a;). (64) 
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It is particularly interesting to analyze this result in the limit of low driving frequency 9q. To 
this end we expand the previous equation up to second order in this parameter. In the ensuing 
result we find that in this limit the mean power contains two different kinds of components 

P l = PT + Pf\ being 

q j. r 2 N r-\-OQ 

PT = -^E^^i) dooRe{ lu (u)D^(u)\, (65) 

c j=l J—oo 1 ' 

if = -Mo^cog^ 5 ) f ^ (w) ^Li^ , ( 6 6) 

where the first component is oc an d it is the leading contribution for low driving frequency. 
Since sin(5; — Sj) is an antisymmetric function, this term satisfies 

JV 

E^ = 0, (67) 
1=1 

indicating that some of the forces deliver, while others receive a net amount of energy. This is 
the essential ingredient for this setup to have an operational regime typical of an engine. In the 



present case, it is particularly appealing The fact that the exchanged energy is of mechanical 
nature. The second component, which is oc 0(0,q), is non- vanishing and accounts for the 
conservation of the energy expressed in Eq. (|15p . This component, thus, describes the net 
amount of energy that is dissipated into the reservoirs in the form of heat. It is important to 
notice that for a driving like (|57p . the exchange component is non- vanishing only in the case 
of a driving with a phase-lags Si — Sj 7^ nn. Another possibility to get finite values for is 
to consider local forces with different amplitudes. It can be in general shown that a multiple- 
parametric driving is a necessary condition to have the exchange mechanism in the present 
setup. 

6. Summary and conclussions 

We have considered a simple model to study heat transport in nanomechanical systems. It 
consists in the usual model for acoustic phonons coupled to phononic reservoirs and generalized 
to include the effect of a time-dependent perturbation in the form of elastic forces acting at 
different places of the structure. We have reviewed a theoretical treatment based on non- 
equilibrium Green's functions, which is analogous to the one previously proposed for electron 
systems under ac-driving. We have derived the equations for the heat current along different 
regions of the setup and shown how can be evaluated in terms of the Green's functions. We 
have also derived expressions to evaluate the power developed by the external ac forces. In the 
limit of weak driving amplitudes and low frequencies, we have analytically shown that the mean 
power done by each force, contains a component that represents the heat dissipated into the 
reservoirs. For a driving with multiple parameters, corresponding for instance to driving forces 
oscillating with a phase lag, there is an additional component, which corresponds to energy 
that can be exchanged between the different forces. This component dominates for low enough 
driving frequencies and enables an operational regime of the setup as a quantum engine. A 
similar mechanism has been identified in driven electron systems like quantum pumps [IS] and 
in coupled quantum capacitors [SO]. The present case has the appealing feature of involving 
mechanical work. 

The present theoretical approach sets a general framework for the study of heat transport by 
vibrational degrees of freedom in driven systems. It can be adapted, generalized and improved 
in many directions, to analyze several interesting situations. In particular, we have considered 
phononic reservoirs, but the present treatment can be easily adapted to deal with photonic 
reservoirs. We have considered coupled atoms, but this treatment can be easily generalized to 
consider coupled molecules, which have internal vibrational degrees of freedom. Improvements to 
model realistic setups with many vibrational modes can also be straightforwardly implemented. 
Extensions of the present treatment to consider the anharmonic effects of Ref. [20] with ac 
driving are also possible. To finalize, it is interesting to mention that a family of Id interacting 
models of cold atoms can be reduced to a Luttinger Hamiltonian, which is basically a model 
of harmonic oscillators. Dynamical evolution of such systems in quenching protocols including 
coupling to reservoirs have been recently considered [38]. The present scheme could be useful to 
study such systems under ac driving. 
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Appendix A. Self-energy corresponding to the coupling to a semi-infinite chain of 
acoustic phonons 

The explicit evaluation of the sum over the normal modes entering the expression of the self 
energy ([29]) can be expressed in terms of a spectral density F a (u) as follows 



being E a>n = \Jk a (l — cos q®)/m a . Notice that the Fourier transform of the self energy is 
perfectly define because is an equilibrium function. 
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